%% dx = (0.125,0.125)
solution = [
1 1.13315 1.28403 1.45499 1.64872 1.86825 2.117 2.39888 2.71828 
1.13278 1.28347 1.45427 1.64783 1.86703 2.11545 2.39732 2.71694 3.07922 
1.2807 1.45097 1.64403 1.86292 0 2.39069 2.70959 3.07132 3.48129 
1.44235 1.63402 1.85143 0 0 0 3.05069 3.45866 3.92071 
1.61515 1.82963 0 0 0 0 3.41547 3.87298 4.39042 
1.79517 2.03389 2.30461 0 0 0 3.79907 4.30552 4.87977 
1.97712 2.24023 2.53852 2.87675 3.25983 3.69363 4.18504 4.74245 5.37436 
2.15447 2.44129 2.76635 3.13475 3.55214 4.02501 4.56083 5.16816 5.85645 
2.31978 2.62865 2.97865 3.37526 3.82467 4.33391 4.91097 5.56486 6.30581 
];
exact = [
1 1.13315 1.28403 1.45499 1.64872 1.86825 2.117 2.39888 2.71828 
1.13278 1.28361 1.45452 1.64819 1.86764 2.11631 2.3981 2.7174 3.07922 
1.2807 1.45122 1.64445 1.8634 0 2.39266 2.71123 3.07223 3.48129 
1.44235 1.63439 1.85201 0 0 0 3.05345 3.46001 3.92071 
1.61515 1.8302 0 0 0 0 3.41926 3.87453 4.39042 
1.79517 2.03419 2.30504 0 0 0 3.80037 4.30638 4.87977 
1.97712 2.24036 2.53867 2.87669 3.25971 3.69374 4.18555 4.74285 5.37436 
2.15447 2.44133 2.76639 3.13473 3.55212 4.02507 4.56101 5.1683 5.85645 
2.31978 2.62865 2.97865 3.37526 3.82467 4.33391 4.91097 5.56486 6.30581 
];
err = [
0 0 -1.33227e-15 -3.55271e-15 -1.33227e-15 7.90479e-14 -9.9476e-14 -5.10703e-14 0 
1.11022e-15 -0.000141557 -0.000246902 -0.000356616 -0.000611128 -0.000864012 -0.000770733 -0.000454413 0 
2.35367e-14 -0.000253238 -0.000414548 -0.000483578 0 -0.00196523 -0.00164104 -0.00090701 0 
-3.55271e-15 -0.000373326 -0.000579833 0 0 0 -0.00276516 -0.00135577 0 
-8.88178e-16 -0.000565926 0 0 0 0 -0.00379722 -0.00155128 0 
5.10703e-15 -0.000300152 -0.000424436 0 0 0 -0.00129849 -0.000853517 0 
2.22045e-16 -0.000132467 -0.000145197 6.5837e-05 0.000117764 -0.00010306 -0.000514167 -0.000399643 0 
5.4623e-14 -4.10007e-05 -4.04089e-05 1.45075e-05 2.12081e-05 -6.03491e-05 -0.000174174 -0.000138758 0 
0 0 -4.61853e-14 1.72751e-13 1.90958e-13 -5.15143e-14 3.64153e-14 0 0 
];
